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Abstract 

The process 7*7 — > u is investigated in the framework of the SU{2) x SU{2) chiral NJL 
model. The form factor of the process is derived for arbitrary virtuahty of 7* in the Euchdean 
kinematic domain. The asymptotic behaviour of this form factor resembles the asymptotic 
behaviour of the 7*7 it form factor. 

A recent analysis of the experimental data on the pseudoscalar mesons production in 
the transition process 7*7 ^ P Q (where 7*, 7 are photons with the large and small 
virtuality, respectively, and P is a pseudoscalar meson) revealed satisfactory agreement with 
the pertrubative QCD (pQCD) predictions [||, ^, ||, § for the asymptotic behaviour of 
the 7*7 — > P form factor at large negative virtuality of one of the photons. However, for a 
correct description of the asymptotic behaviour of the transition form factor it is not enough 
to apply only the pQCD technique, and it is necessary to use the nonpertrubative sector 
for determining a distribution amplitude (DA) of quarks in a meson. In the lowest order of 
pQCD, the light-cone Operator Product Expansion (OPE) predicts the high behaviour 
of the form factor as follows ||2|, |^ : 

Fpid, 4,P' = 0)|q._ = Jp M ^ + o (^) + o (^) , (1) 

with the asymptotic coefficient given by 

1 



where fp{x) is the leading- twist meson light-cone DA normalized by 

/ dx ipj,{x) = 1; (3) 
Jo 

as is the strong coupling constant; and ui are, respectively, the total virtuality of the 
photons and the asymmetry in their distribution 

Q"" = -{ql + ql) > 0, and u;={ql-ql)/iqf+qj), \u;\^l, (4) 

where qi,q2 are the momenta of photons; fp is the meson weak decay constant (for the pion 
= 93 MeV). 
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Unfortunately, the determination of ^p{x) is a rather nontrivial problem and can not be 
performed in the framework of only pQCD. At asymptotically high Q^, DA is ipp''^-^'^^^(x) 
= 6x(l — x) and j^^^p^ (||^| = i) = 2. The fit of the CLEO data||l| for the pion corresponds 
to J^^^'^ (|w| ~ 1) = 1.6 lb 0.3, indicating that already at moderately high momenta (Q^ = 8 
GeV^) this value is not too far from its asymptotic limit. 

In our previous work [^, it has been shown that the calculation of a process 7*7 —> P 
amplitude in the framework of a chiral quark model of the NJL type is free from difficulties 
connected with the necessity of determining DA. Our results for the form factor asymptotics 
agreed with experiment Q, and with the predictions made in QCD sum rules @, I, g and 
the instanton induced quark model (IQM) [^. Also, both the NJL and IQM model allow 
one to derive the meson DA. 

Here, we continue the investigation started in Q and consider the scalar isoscalar meson 
production through the process 7*7 a where a is associated with the lightest scalar 
isoscalar state /o(400 — 1200) [l^, and the photon 7* is off-shell. Direct observation of a is 
hardly possible; however, in the low-mass region it can show itself as a resonance in the pion 



pair production [11|. 



In our paper, all the calculations are performed in the framework of the SU{2) x SU (2) 
NJL model. We derive a gauge- invariant expression for the amplitude of 7*7 — > a and 
determine the form factor of the process for which the asymptotic behaviour is studied at 
large virtualities of 7*. Notice that the asymptotic behaviour of the 7*7 — > a and 7*7 — > vr 
form factors is similar. 

The interaction of mesons with u and d quarks is described by the following quark-meson 
Lagrangian: 

L{q, q, A, a, tt) = q{x){id - M - eQA{x) + g{a{x) + ij5T7r{x)))q{x), (5) 
where q and q are the u and d quark fields 

"f^) - ( ) ■ (^) 

M is the constituent mass matrix, M = diag(mu,m(i) (m„ ~ rrid w m); A = A^7^ is the 
photon field; a stands for the scalar isoscalar meson; tt is the pion triplet, r are the Pauli 
matrices ; e is e.m. charge e^/47r = 1/137; Q = 1/6(1 + St^) is the charge operator; and g is 
the constant of a quark-meson interaction defined by the relation g~^ = 4I2 where 

3 /■ 4fc d{A^ - e) 

^ (27r)4 J (A;2 + m2)2 ' ^' > 

The subscript E in d'^k means that the integration is performed for momenta in the Euclidean 
metric. The cut-off A eliminates the UV divergence. 

Using the Goldberger-Treiman relation g = m/ f-,^ and the relation gp = ^/^g [^, where 
gp Q.l is the constant describing the decay p — > 27r, we determine the constituent quark 
mass m = 234 MeV and A = 1050 MeV. 

Let us note that the 7*7 a amplitude also describes a two photon decay of a. The 
amplitudes for radiative decays of scalar mesons (5) were obtained in The amplitude 
of the process — > 77 has the form 

A = e^m^^e,{q,)eM2\ (8) 
where e are polarization vectors; qi and q2 are photon momenta. 



^Following iQl, we do not take into account u — ai transitions. Therefore, = 9a 
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Up to the one-loop, the process is described by the diagrams in Fig. |l] that give for the 



tensor T^'^^^ 



where 



-^5^77 ~ 5^77 "I" -'5^77' V^-* 



ri^;; = Cs{qU'^-g^''{qi-q2))Ju.{quq2), (10) 
tS;; = CsB^-{q^,q2), (11) 
and the transversahty condition efj_{qi)q^ = £1/(92)92 = is taken into account. Here 



(27r)4 J D{Q)D{q^)D{-q2) 
D{q) = m^-{k + qf. (14) 

The constant Cs is to be specified for each type of scalar mesons. In the case of a, it is 

40 

Ca = (15) 

The term T^^^-y has a gauge-invariant form. The Pauli-Villars regularization is necessary 
to restore the gauge invariance of Tj^}^^^ which leads to the following form of B'^'^ (see 
Appendix for details) for 9i = 9I = 0: 

B'^^((/l,92) = 2(j^^-g^'''\m\jrr,{<d,{))-J„Mi:q2))- (16) 



, {qi ■ q2 

In our paper, we apply the Pauli-Villars regularization with one subtraction for an arbi- 
trary function Y{m?') as follows 

yp^(m2) = y(m2) - y(M2) (17) 

Formally, as all integrals we calculate here converge, the regularization can be lifted off by 
reaching the — > 00 limit. In our particular case, the regularized expression will, however, 
differ from the non-regularized one, the constant term, which violates gauge invariance, is 
thereby cancelled. 

In the case of a production through 7*7 — > a, the amplitude can also be divided in two 
parts 

rpfiV _ rp(l)lJ,l' ,rp(2)lJ,U 

-^7*7 — s-cr -^'y^'y — ^cr ' -^7*7 — >cr^ \ "^J 

where the first term is again gauge-invariant 

T^l\Z. = CAqU'^-g^-{qi-q2))Jm{qi,q2), (19) 

whereas the second one becomes gauge invariant after the Pauli-Villars regularization is 
implemented and then lifted off. Moreover, as one of the photons is off-shell, B^^ acquires 
an additional term proportional to 9^ (see Appendix): 



(20) 



\[qi-q2) 



2m^(J„(0,0)- J^(9i,92)) 

■ q2) / 

^ f^4, W 1 1 



2(91 • 92) (2^)4 J D{q{) \D{-q2) D{0) 



After the gauge invariance of the amphtude is restored, we define the process form factor 
^(91,92): 



where 



rplJ,U 



40/^5' 



(211 



2{q, . (2vr)4 



Jm{qi,q2) 



2m^ 



{qi ■ 92) 
1 / 1 



{Jm{qi,q2) - Jm(o,o))- 
1 



Diqi) \Di-q2) DiO) 



(22) 



After the change of variable k = k' — {qi — (?2)/2 (in the fohowing we omit the prime symbol), 
the form factor takes the form 



F{qi,q2) 



2m^ 



5/^ 
67r2 

1 



{qi ■ q2) D{Of 



d'^k 
+ 



1 - 2mV(gi • 92) 



D{-p/2)D{p/2)D{-{q^-q2)/2) 
ql 1/1 



2(gi • q2? D{p/2) \D{-p/2) D{-{q, - q2)/2) 



where p = qi+q2- It is convenient to consider F(qi, 52) in terms of p'^ (notice that p^ = M^) , 
qf, and Further, we restrict ourselves to the case of ^2 = 0; thus, we have two independent 
variables and qf. The kinematics of the process under investigation corresponds to large 
negative qf, and we introduce, for convenience, a positive magnitude defined as = —qf- 
After that, the form factor F{qi, ^2) can be considered as depending on p^ and 



F{p\Q'')=F{qi,q2) 



Ig|=0;g2<0> 



(23) 



where 



F{p\Q^ 



67r2 



4m2 \ , 2 ^2n 



2Q2 



(p2 + Q2)2 



{MAip')-NAip^Q')) 



Kk{p\Q^) 
La 
Ma(/) 

Na{p^Q') 



^d^k 



1 



m2 D{-p/2)D{p/2)D{-{q, - q2)/2) ' 
^d^^ 1 



i7r2 Z)(0)3' 
^d^^ 1 



i7r2 D{-p/2)D{p/2) ' 

^d^^ 1 

m2 Z)(-(gi-(72)/2)I?(i^/2)' 



(24) 

(25) 
(26) 
(27) 
(28) 



In all integrals, we implement the UV cut-off A that constrains quark momenta to the 
domain where chiral symmetry is spontaneously broken and the bosonization of quarks takes 
place (see (|^). 

An analogous method has been used in 0]. One can compare ( p4|) with the pion form 
factor obtained in for / 0, 7^ 



J ZTT^ U 



{-p/2)D{p/2)D{-{qi-q2)/2)' 



(29) 
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where is a constant. In 0, the fohowing expression was obtained for Kx{p'^,Q'^) after 
introducing Feynman parameters, integrating over the 
angles, and changing the variables /c^ = u: 



duu 



m? + u — ^ 



(30) 



1 



X dx 







+ 



Here 



b = m'^ + U + ^xQ^ - ^ (1 - 2a;) 



a± = 2uxQ^ {x±{l- x)) - (1 - 2x) np^ 
Analogously, one can write the integrals N\{p'^ ,Q'^),M\{p'^) as follows: 



Na{p^Q^) 
Ma{p^) 



A2 

2 duu dx 





A2 



duu 



{ bo ibo + y/bl- a+o 



Here 



6o = 6U'^o = m +u - -p , a+o = a+U-*o = -wp 



(31) 

(32) 
(33) 

(34) 



Now we can calculate the asymptotics for F{p'^,Q'^) when oo. Using the approx- 

imation described in Q and expanding in a series of l/Q'^, we obtain for K/^{p'^,Q'^), La , 
MA(p2) and Na{p'^,Q'^) at = q 



A2/m2 



-^A = — 5- 



Ma(p' = 0) 



du u 
(l+n)3' 



AV»^ 



A^a(p' = 0,Q2) 



2m^ 



duu 
(l + u)^' 

AVm2 



y dnln(l+n)+0(^ij) 



(35) 
(36) 
(37) 
(38) 



Note that the last term is of the order of As a result, the first three terms give us the 

following asymptotics for the form factor F(p'^ = 0, Q^): 



F(p2 = o,Q2) 
Js 



Q 



1 



3.28 + 0.96 - 2.22 = 2.02. 



(39) 
(40) 
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We also calculate numerically Q'^FijP', Q"^) for intermediate values of and for a different 
choice of p^. The curves drawn for = 0, = = (0.45 GeV)'^ are shown in Fig. |2|. The 
value M(j = 450 MeV is consistent with recent experimental [|l3| and theoretical data. A 
comparison of the obtained results for the function Q'^F{p'^,Q^) with an analogous function 
for the pion [0] shows similarity in their asymptotic behaviour. However, some differences 
take place at intermediate values of Q^. Indeed, the pion function grows monotonically in 
the whole region of Q^. At the beginning it grows rapidly and after > 2GeV^ it slowly 
approaches the asymptotic value. In the case of the a meson the function Q"^ F {p^ , Q"^) also 
experiences a rapidly growth up to approximately 2GeV^; after that it is slowly decreasing 
to an asymptotic value. Asymptotic values for these functions are practically equal to each 
other. 

In this work, we have shown that in the framework of the SU (2) x SU (2) chiral quark 
model of the NJL type, it is possible to describe the behaviour of the 7*7 a form factor in 
a wide region of Q^. The exact expression for the form factor of the process is obtained, and 
its asymptotic behaviour is investigated. A comparison of the vr (see 0) and a form factors 
reveals similarity in their asymptotic behaviour, which may be understood as a consequence 
of chiral symmetry. 

Our results can be useful in investigations of the processes where pion pairs in the S'-wave 
are produced in two photon collisions [15|. Indeed, in these processes, as a rule, the a-pole 
diagram plays the dominant role Some data on the pair pion production is already 

available Also, one can find a discussion on this topic in [p^ j. 

Notice also that the information concerning the amplitudes 7*7 — > 7r,7*7 — > a can allow 
us to calculate corrections to the muon anomalous magnetic moment from the processes 
7*7 — > vr ^ 7*7*)7*7 — > o" — > 7*7*) where 7* interact with the muon. Last year, this topic 
has been discussed in various papers (see 

Further, we plan to calculate the transition form factor for two off-shell photons with 
arbitrary virtualities in the framework of both the NJL and IQM model, and it allow us also 
to define DA. 

The authors thank A.E. Dorokhov, S.B. Gerasimov and E.A. Kuraev for fruitful discus- 
sions. The work has been supported by RFBR Grant 02-02-16194. 
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Appendix 

Let us rewrite the expression for JmiQi,Q2) (P^, using the Feynman parameterization 

1 1-x 

2i f , f , f d!^k 



dx / dy 



{2iiY J J J {m"^ — k"^ — 2k{qix — q2y) — q^x)^ 


Shifting the variable k ^ k — {qix — q2y) and integrating over momenta and y, we obtain 

1 

^ , , 1 1 f dx ^ ( m? — p^xil — x)\ , 

Jm{qi.q2)=-j^ir, r / —In ^ \ ( . 42 

(47r)^ 2(gi • ^2) 7 X - q{x[\ - X) ) 





Let us rewrite the expression for B^^ (see (|13|)), using the Feynman parameters 



dx / dy d k 



(27r)^ J J J {m? — k"^ — 2k{qix — q2y) — qfx)^ 


Again, after the change of variables, we have for B'^'^ 

1 l-x 

2i f f f A^^ 
B^^qu,.) = -j-y,jdx j dyjd-k _ _ _ _ ^^^^ ^^^^^^3 , 



A^''' = g^^^irr? -k^ + 2k- {qix - q2y) - qW + 2(^1 • q2)xy) + 

+4(A;^r - k^'iqix - q2yY " {qi^ - q2yTk'' + (44) 
+ {qix - q2yT{qix - q2yY)- 

The integration over momenta gives 



1 l-x 

gflu r r ^2 _ ^2^2 _^ 2{qi ■ q2)xy 



B^'^{q,,q2) = ^,jdx I dy 



' w? - qix{l — x) - 2{qi ■ q2)xy 



#1 ]dx Jdy , 2 ^ "\ ^ ^- (45) 

(4vr)^ J J - qlx[l - x) - 2{qi ■ q2)xy 





Let us denote the denominator in ( p5| ) as Z 

Z = — q\x{l — x) — 2{qi ■ q2)xy. (46) 
Using this equation, we substitute the mass squared in the numerator of (^) 

1 l—x 

9^" / , /■ , ^ + 4(91 • q2)xy + qlx{l - 2x) 



B^''{qi,q2) = ^ j dx j dy 



Z 
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Now it is easy to separate the gauge-invariant and noninvariant parts in 

1 l-x 



{qi,q2) = / dx / dy— (48) 





1 l-x 

9^" (' [ , Z + qlx(l - 2x) 





The first term in ( ^8|) is gauge- invariant, while the other is noninvariant and can also be 
divided in two terms 

1 l-x 1 l-x 

^ Ux [ dy+-^ I dx f dy^^^^^. (49) 



(47r)2 J -^"^ J ' (47r)2 , 





The first term in ( |49| ) is a constant and is cancelled by the Pauli-Villars regularization. 
Integrating over y in the second term, we obtain 

1 

'm? - {qj + 2{qi ■ q2))x{l - x) " 



(50) 



(47r)2 2{qi • (72) 7 V — qfx{l — x) 



After the substitution 1 — 2x = z, one can see that this part equals zero 



(47r)2 4(gi • ^2) 7 V 4m2 + g2(i_22) 



-1 

Therefore, the expression for B^^'^ becomes gauge invariant after the Pauli-Villars regulariza- 
tion is implemented and then lifted off 

B>"'{qi,q2) = iqU^-{qi-q2)9ninm')-HM'))\M^^oo, (52) 

1 l-x 

1 f f Axy 



I{m ) = J J ^^~Z~' ^^^^ 



Integrating over y in I{m ) we obtain 





2^ 



1 



^, 1 1 1 m? f dx , f m? — p^xil — x) 



(47r)2 (gi • 52) (47r)2 (gi ■ 92)^ 7 x \w? - qlx{l - x) 



1 

''^ 'd2;(l-x)ln ^ — ^-2-7- f . (54) 



(47r)2 (qi ■ (72)^ J \?Ti2 — qfx{l — x) 





.2\ 



Now, let us show that I{m ) can be expressed through the integral 

' /■d4^. ^(0)-^(-92) 



(27r)4 J D{0)D{qi)D{-q2) 

KxTdyfd'k (^-^^ 



(27r)^ J J J {nn? — kp' — 2k{qix — q2y) — qfx)^ 


j:^/rfxin( '";-''^;;-"i ). (55) 

(47r)2 7 - qfx{l - x) J 
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Noting that the fohowing identity is satisfied: 

1 1 

1 f ( m'^ -p'^x{l- x) \ 1 /■ , , ( - x) \ 

dx[l - xjln — 5 1^— = , ,„ / dxln — ^ , (56) 



(47r)2 y \m? — q\x{l — x) J 2{A'kY j — q\x{l — x 



we obtain for (^5|) 

^ (' D{0)-Di-q2) _ 2 } , f m^-p^x{l-x) \ 

(2vr)4y " D{0)D{q,)D{-q2) ~ (4^)2 7 ^"^^^ W - qfxil - x) J ' 



Using (42) and (|57|), we rewrite I{m?) as follows: 

2\ 11 2m? 

' ^^^fe.^^°^r^i7'^^- (58) 



2(gi • 92)2 (27r)4 7 L>(0)Z)(gi)D(-g2 
Formally, the first term in ( |5^ ) can be rewritten as 

2m^ 



because 



Jm(0,0) (59) 



(fk _ 1 
(27r)4 y D(0)3 ~ 2(4^2j„2 



•^-(0' 0) = -TKZ^A I 7^ = nr. .,2 2 " (60) 



Finally, the regularized expression for B'^'^{qi,q2) has the form 

/ 277)2 

i?^"(9l,'?2) = («-(5l-'Z2)5^'^) -7 Mqi,q2) + -, rJ^(0,0)- 



(91 • 92) {qi ■ 92) 

11 ^ [ D{0)-D{-q2) 



2{q, . q^Y (27r)4 j ^ ^ D{Q)D{q^)D{-q2) ) ' ^^^^ 

Equation (^) is expressed in terms of formal integrals. This gives us an advantage of further 
implementing a regularization different from the Pauli-Villars scheme. 
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?2, 



92. f,^ 



Figure 1: Diagrams contributing to the amplitude of the process 7*7 cr 




Figure 2: The form factor of the process 7*7 a muhiphed by at from to 10 GeV^. 
The thick curve is for = (0.45 GeV)^ and the thin is for p"^ = 0. We compare these curves with 
the theoretical limit for the pion transition form factor: 2/^^ 
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